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Abstract 

We propose a dual thermodynamic description of a classical instability of generalised 
black hole spacetimes. From a thermodynamic perspective, the instability is due to 
negative compressibility in regions where the Casimir pressure is large. The argument 
indicates how the correspondence between thermodynamic and classical instability for 
horizons may be extended to cases without translational invariance. 



1 Introduction 



The intriguing connection between classical and local thermodynamic instability [1, 2, 3, 4] 
has partially been behind a renewed interest in the Gregory-Lanamme instability of black 
branes [5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]. The connection states that translationally 
invariant horizons are classically stable if and only if they are locally thermodynamically 
stable [1, 2, 3]. The objective of this work is to suggest a dual thermodynamic description of 
a recently discovered black hole instability The novel aspects of the connection presented 
here are that the horizon is inhomogeneous and that the thermodynamic instability is due 
to a Casimir effect. 

Local thermodynamic stability requires the Helmholtz free energy of the system to be a 
concave function of temperature and a convex function of volume. The known instabilities 
of black branes are due to a negative heat capacity, Cy = —Td 2 F/dT 2 \y < 0. In contrast, 
the instability in this paper will be due to a negative isothermal compressibility K^ 1 = 
Vd 2 F/dV 2 \ T < 0. 

We work within the context of the AdS/CFT correspondence [16, 17, 18]. Amongst 
other things, the AdS/CFT correspondence relates the physics of black holes in AdS^+2 
to a thermal field theory living on the 'boundary' S l x S d [18, 19]. In the conclusion we 
note that it should be possible to extend the ideas presented here to general inhomogeneous 
horizons, not necessarily embedded in an AdS-like geometry. 

Section 2 reviews an instability of generalised black holes [20, 21] with a negative cos- 
mological constant [22]. We further recall how the criterion for classical instability of the 
black hole translates into dual field theory language. 

Section 3 calculates the energy momentum tensor of the thermal field theory dual to the 
black hole spacetime using holographic renormalisation [23, 24, 25, 26, 27]. The thermal field 
theory has an ultrastatic background S 1 x M, with M. an arbitrary dimension d Einstein 
manifold. There exists an ambiguity in the energy momentum tensor due to local curvature 
polynomials of mass dimension (d + l)/2. Generically, these produce an inhomogeneous 
Casimir energy and pressure. We show how these terms may cause a thermodynamic 
instability. This thermodynamic instability appears to be a suitable dual description to the 
aforementioned classical instability of the black hole bulk spacetime. 

Section 4 illustrates the previous considerations in the case when the horizon of the 
black hole is given by a Bohm metric on S 5 [28, 21]. The thermal field theory in this case 
has background S 1 x S^ 6hm . 

The last section is the discussion. The argument we have presented is not rigorous and 
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we clarify what would be needed for a watertight discussion. Further applications of the 
Casimir effect to horizon instabilities are suggested. 

2 Instability of generalised AdS black holes 

We start by reviewing the stability of generalised black holes [20, 22] with a negative cos- 
mological constant. Generalised black holes take the form 

ds 2 = _ f{r)dt 2 + ^_ +r 2 ds 2^ (1) 

where f(r) = 1 — (1 + r+/L 2 )(r + /r) d ~ 1 + r 2 /L 2 . The event horizon is at r = r + . If ds 2 d 
were the round metric on S d , then the spacetime would be the usual Schwarzschild black 
hole in AdS space. However, the d + 2 dimensional vacuum Einstein equations allow the 
horizon metric ds 2 , to be any d dimensional Einstein metric on a Riemannian manifold A4 

Rab = {d~ l)9ab ■ (2) 

The resulting spacetime is called a generalised black hole. The asymptotic geometry of the 
background has now changed and the dual thermal field theory lives on S* 1 x M. 

The stability of generalised black holes has recently been investigated [22, 20, 21]. For 
large, r + /I>l, generalised AdS black holes one finds 

A L <~xO(l) ^ instability, (3) 

Li 

where is the minimum eigenvalue of the Lichnerowicz operator on M. Recall that the 
Lichnerowicz operator acts on symmetric rank two tensors as 

(A L h) ab = 2R c abd h d c + R ca h c b + R cb h c a - V c V c h ab . (4) 

Note from (3) that the minimum Lichnerowicz eigenvalue needs to be very negative in order 
for instability to occur. Examples of Einstein manifolds with large negative Lichnerowicz 
eigenvalues are the Bohm metrics on S 5 . . . S 9 and on products of spheres [28, 21]. We will 
use these metrics as examples below. 

The criterion may be translated into field theory language using the standard AdS/CFT 
dictionary to give [22] 

T 2 < -X L x 0(1) instability, (5) 

where T is the temperature of the field theory. Thus, the AdS/CFT correspondence predicts 
a critical temperature in the dual field theory, T 2 , ~| Xl \. 
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The Lichnerowicz operator is not a particularly natural object to consider in field theories 
where gravity is nondynamical. On Einstein manifolds with positive curvature the minimum 
Lichnerowicz eigenvalue is related to the Weyl curvature. Roughly, one expects [21] 

a| ~ c a bcdC abcd \ max . ~ yi 2 , (6) 

where we use 9^ to denote a typical magnitude of curvature of the geometry M., and 
C a bcdC abcd \max. is the maximum value taken by C a \ ia iC ahcd on the manifold. Putting these 
statements together, the critical temperature becomes 

t£~<r. (7) 

Note that the Ricci scalar, R, is fixed, so in the limit of large negative Al we have 

This hierarchy is important for consistency because the large black hole limit r + /L » 1 

translates into the field theory condition T » R. 

In the next section we present a potential field theory instability occurring at this critical 
temperature. 



3 Casimir instability of field theory thermodynamics 
3.1 Holographic energy momentum tensor 

The thermodynamics of the dual field theory in the strong coupling regime may be studied 
using holographic renormalisation [23, 24, 25, 26, 27]. This is the regime dual to the weakly 
curved gravitational description in the bulk and is therefore the regime in which we expect 
to find a field theory instability corresponding to the black hole instability. 

For concreteness we will work now with seven dimensional black holes and a six dimen- 
sional boundary. This allows the problem to be embedded into the AdSj x S A version of 
the AdS/CFT correspondence and will allow us in the next section to use five dimensional 
Bohm metrics as concrete examples. We will work from now on with a Euclidean signature, 
to describe finite temperature physics. 

The expectation value of the energy momentum tensor of the dual theory is given by 

2 dS 



< y/w > = 



(8) 

rcnormaliscd 



\/det g dg^ 

where gw is the metric on the boundary of spacetime. The gravitational action, S, has 
three parts: the Einstein-Hilbert action, the Gibbons-Hawking boundary term and the 
counterterm action. The counterterm action is defined on the boundary of spacetime but 
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differs crucially from the Gibbons-Hawking term in that it is a local functional of intrinsic 
boundary curvatures only. This implies that it does not affect the bulk equations of motion. 
The counterterms are added to cancel divergences in the action due to the diverging volume 
of the spacetime at infinity. 

The calculation of the dual energy momentum tensor for a general spacetime with nega- 
tive cosmological constant and a six dimensional boundary has been performed, for example, 
in [23] . It is straightforward to apply their formulae to the generalised black hole spacetime 
(1). Write the boundary metric as 



g^ v dx^dx u = dr 2 + g a b dx a dx b 



(9) 



We have fx, v running over all the boundary coordinates and a,b . . . running over the spatial 
directions of the boundary. We find 



< Too > = 
<T ab > = 



3L 5 



87rG 7 rg 

3L 5 
8vrG 7 rg 



6T 4 



L 2 



48 



lr 



1 



goo + Tbolf.c. , 

9ab + |f c j 



(10) 



L 2 J 48_ 

where f.c. stands for 'finite counterterms', described below, and ro is an arbitrary scale 
introduced during the renormalisation. Dependence on vq will disappear shortly. The 
result (10) may be reexpressed in terms of field theory quantities by relating the radius of a 
large black hole, r + /L S> 1 , to the field theory temperature, using the AdSi x S A relation 
between the AdS length and Newton's constant and introducing the Ricci scalar of g a b\ 



This gives 



<T 00 > 



<T ab > 



2tt roL 

N 3 
3^3 



1 

G~7 

f y T6 



3^U 



R — — o . 

r 2 



(11) 



5-2 3 



R d 



9oo + Too |f 



c. ' 



2vr 



5-2 3 



R< 



9ab + T a b If .< 



(12) 



The Ricci scalar term is in fact negligible because the large black hole limit implies T 2 3> R. 

The origin of the finite counterterm contribution is as follows. The counterterm bound- 
ary action consists of curvature scalars. The terms in the action that are scalars of mass 
dimension less than six diverge as the boundary is taken to infinity and are used to cancel 
the divergences of the bulk gravitational action. Scalars of mass dimension greater than 
six are negligible as the boundary goes to infinity. Scalars of mass dimension precisely six 
remain finite. Thus, one is free to add these scalars to the countertem action. 
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The same ambiguity also exists in the dual field theory. The coefficient of a mass dimen- 
sion six curvature scalar in a six dimensional action is dimensionless and should therefore be 
viewed as a renormalisable coupling. The value of the coupling is thus generically not zero. 
Its value is not determined by the theory but needs to be observed 'experimentally'. The 
contribution of these terms to the action might be thought of as a position dependent cos- 
mological constant. The variation of these terms with respect to the boundary metric gives 
the contribution to the energy momentum tensor that we have included above in (10) and 
(12). Due to the fact that there this constitutes a temperature independent contribution to 
the energy and pressure, we call the presence of these terms a Casimir effect. 

3.2 Thermodynamic instabilities - general argument 

Before considering explicitly the possible form of the finite counterterm contribution to 
the energy momentum tensor, we will see how such terms can produce thermodynamic 
instabilities. The following section offers a concrete example. Recall that we are interested 
in manifolds with a curvature scale much larger than the constant Ricci scalar » R. The 
picture one has in mind is the following. 




Figure 1: Manifold with two curvature scales, Riem | 

max ^ 

9\ 2 > R 2 . 

In the region where the curvature scale 9\ is large, V\ » R, the pressure in (12) will be 

P ~ N 3 [T 6 ± 9* 3 ] . (13) 

We are neglecting positive order one coefficients. We have taken the N 3 term as an overall 
factor because one would expect the Casimir contribution to be proportional to the number 
of degrees of freedom. 

If the temperature is sufficiently high, T 2 3> £H, then the field theory thermodynamics is 
just that of black body radiation of 0(N 3 ) species. However, if T 2 < 9\, then the Casimir 
pressure is dominant in the regions of large curvature. This is consistent with the large 
temperature approximation we are making, T 2 » R, because of the hierarchy of curvatures 
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5H > i?. Thus, the qualitative thermodynamic behaviour in the strong curvature regions 
appears to change at around Tq ~ 9\. 

The curvature 91 has an associated lengthscale Ly^ = 9i -1 / 2 . This lengthscale gives the 
size of the strong curvature region. The lengthscale associated with the wavelength of ther- 
mal excitations is = T -1 . In the lower temperature regime T 2 < T% ~ 5H, the thermal 
wavelength is larger than the size of the strong curvature region: > L<r. Therefore the 
bulk of the thermal spectrum is not excited in this region. Effectively, the region is not 
in thermodynamic equilibrium with the rest of the radiation. To first approximation, the 
strongly curved region is at zero temperature. In this case, let us consider the behaviour of 
the strong curvature region as an independent system. The volume of this region is 

/ 1 \ V2 

V~l&= (<= • (14) 



Writing the pressure (13) in terms of the volume of this region 

P ~ ± (v) • (15) 

Consider the isothermal compressibility 
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The compressibility may also be expressed in terms of the Hclmholtz free energy by using 
the Maxwell relation P = -dF/dV\ T 

= " PI" (17) 

Thus where the compressibility is negative, the free energy is not a convex function of the 
volume and the system is thermodynamically unstable. 

We see that there will be a thermodynamic instability at T 2 = T 2 , ~ $H if the sign of the 
Casimir pressure in (13) is negative. We will see below that in practice the contributions to 
the pressure come with both positive and negative signs. The present argument indicates 
the generic possibility of instability, rather than a necessary instability. The localisation 
of the instability to the strong curvature regions is mirrored by the classical black hole 
instability, where the unstable mode is concentrated in these regions [21]. 

The dual operator to the unstable mode [22] is a linear combination of spatial com- 
ponents of the energy momentum tensor. Below the critical temperature one expects the 
vacuum expectation value of this operator, which are components of the pressure, to shift. 
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This suggests that there is a new stable phase where the curvature does not give the domi- 
nant contribution to the pressure. 

This instability should be compared with previous matchings between dynamical and 
thermodynamic instability [1, 2, 3, 4]. The instability in these works is due to a negative 
heat capacity 

— (I) -(SX- 

3.3 Contributions to the Casimir pressure 

We can be more precise about the curvature terms contributing to the energy-momentum 
tensor in the strongly curved regions. In six spacetime dimensions, there are 17 curvature 
scalars that have mass dimension six [29]. These will generically be present in the action 
and their variation will contribute to the energy momentum tensor. However, not all of 
the contributions will be large. The spatial background is Einstein, and consequently has 
constant Ricci scalar, R, and Ricci tensor. Some of the tensors will be 0(R 3 ) or O(R 2< 0\) 
or O(R0\ 2 ) instead of C(9^ 3 ). These will be negligible in the region where 9\ » R. 

Given that the field theory does not specify the coefficients of the curvature terms in 
the action, there is not much to be gained at this stage from starting with the action. 
Instead of calculating the variation of the 17 curvature scalars, we will start by directly 
considering all possible rank two tensors that have mass dimension six. Although there are 
many such tensors [29], if we consider a background of the form S 1 x Ai, with Ai Einstein, 
then there are only eight independent rank two tensors with mass dimension six that are 
0(m 3 ). Explicitly, 

8 

T ^i f .c. = E^ + °(^ 2 )' ( 19 ) 

for some coefficients ki and where the tensors K$, are 1 

^ ^ — ^RP^ V (ju^ y^Rpqrs i ^ fiR^"^ V i/Rpqrs j V S .R^ pX7 Rpqrvi R^ Rpq ^.Rrstv i 

RPI Rp r^iRqtsv i RP^ Rpqr Rs/itu i Vt-R^ V RpqrsQjiv i Rf"^ Rp r Rqtsu9fJ-f} • 

(20) 

All other linearly independent rank two tensors with mass dimension six either vanish or 
contain factors of the Ricci scalar when considered in an Einstein background. This may 
be worked out from the basis of tensors given in [29]. The S 1 factor in the metric does not 
J For a general six dimensional metric there is a ninth independent tensor R pqrs R vq tuRrs tu g^v , but this 
is not independent in five dimensions, or in the case S 1 x M. Thus in our backgrounds it is not C($H 3 ). 
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contribute to any of the curvature tensors and so all the summations involving curvature 
tensors are over the five spatial directions - compare with the notation of equation (9). 

Not all of these terms are obtainable from varying an action. We must impose two 
constraints. Firstly, diffeomorphism invariance of the action implies that the energy mo- 
mentum tensor is transverse (the Noether constraint) because under the diffeomorphism 

"-W^ 7 ^*"/^) 4 * =* v ^r v ^ =0 - (21) 

Secondly, in six dimensions an action given by scalars with mass dimension six is invariant 
under constant Weyl rescalings. This implies that the trace of the corresponding energy 
momentum tensor must be a total derivative. Under Sg^ = ag^: 

I -^-ag^dSl = 5 a S= f V^aJ^dn = - f aV^dfl g^^- = T = -V ^ , 

J og^i; J J ^g^iv 

(22) 

for some vector J^. This total derivative contributes to the non- universal part of the Weyl 
anomaly [30, 31]. 

The trace constraint is the easier to apply. The possible candidates for are rank one 
tensors with mass dimension five or indeed zero. The only relevant nonzero possibility is 
J M = R pqrs V f+Rpgrs [29]. Considering the eight tensors in (20) it follows that there are seven 
admissible combinations 

Ik^+kW K<V+K<V K<U + -K<n K®--K<® K<$ + 4K®\ 

(23) 

where, in taking the trace, any terms arising that are not Offi 3 ) are ignored. It is assumed 
that any contributions to the trace smaller than 0(JK 3 ) may be cancelled using other small 
terms. The footnote of the previous page is important for this calculation. 

Now consider the divergences of these seven tensors. We need combinations of the 
tensors in (23) that have vanishing divergence. The appendix indicates how one goes about 
finding the combinations systematically. The result is four linearly independent tensors 

H {i) _ K (4) _ 2R (5) _ K {6) , l K (8) 

H (2) = _ R (1) + R (2) _ 16R (5) + i R {%) 

H$ = -K$-SK$ + \k$-\k$. (24) 
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Let us summarise what we have just achieved. We wanted to calculate possible contri- 
butions to the energy momentum tensor coming from a mass dimension six curvature scalar 
in the action. Further, we were only interested in terms that are C(9^ 3 ) and hence give 
the dominant contribution in the strongly curved regions. These contributions had to be 
rank two tensors with mass dimension six. We have shown that the only tensors satisfying 
all these conditions are the four in equation (24). We may now calculate these terms in a 
concrete background. 

4 Dimension six tensors of Bohm metrics on S 5 

The Bohm metrics [28] serve as an illustration of the ideas we have introduced. We will 
consider Bohm metrics on S 5 . The metrics have the form 

ds\ = de 2 + a{efdn 2 2 + b{efdh 2 2 , (25) 

where dti 2 , and dQ 2 are round metrics on spheres S 2 and S 2 , with polar coordinates (a, j3) 
and (a, (5) respectively. The coordinate 9 has range to 9f. The Einstein equations imply 
some nonlinear differential equations for a(9) and b(9). The boundary conditions for topol- 
ogy S 5 are that a(0) = 0,5(0) = bo and a(9{) = ao,b(9{) = 0. In fact, the solutions to the 
equations on S 5 have a{6) = b(9f — 9). There is a discrete infinity of solutions, at specific 
values of bo. Solutions exist for arbitrarily small &o- As &o —> the metric develops conical 
singularities at 9 = 0, 9 = 6{. Thus &o gives the resolution of the singularity. We expect to 
find 9t - b Q 2 . 

We can evaluate the four contributions to the energy momentum tensor of (24) at the 
points 9 = and 9 = 9{ using Taylor series expansions for the functions a(9) and b(9) [21]. 
The results for bo <C 1 are contained in Table 1 below. One may obtain the values at 9 = 9{ 
by letting H ( J )a <-> H% )& , HfP and H^ )e remains the same. 

We are working in units where the Ricci scalar R = 20, so the condition bo <C 1 is 
the introduction of a new scale: 6q 2 ~ S> R. The first observation we can make about 
Table 1 is that the contribution to the pressure is Offi 3 ) as expected. The pressure in each 

(i) (i) 

direction is given by Px = kiH x . The second observation is that contributions to the 

(2) 

pressure from Hj^J have different signs for different directions. Thus, whichever the sign of 

(2) 

the coefficient, Jt2, of H^J in the energy momentum tensor, one of these pressures will be 
negative, contributing towards a thermodynamic instability as outlined above. Of course, 
all the pressures could still be made positive by adding a larger positive contribution from 
one of the other tensors. 
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Tensor components 


At = 


±± 6 ±±a ±J - (3 "a " p 


32 
276g 


l^^Ol " ft !"o J 




f 64 128 256 | 
1 276 s , ' 276[j ' 276[j j 




{H®°, H^o = Hff>, Hf* = Hft) 


< 


f 1 32 176 1 
{ U \ h 966 ' 276 6j 


> 


(rrWe o-(4) a _ rr(4)/3 ff (4) fi _ rrW^l 


{ 


64 256 320 
276g ' 276g ' 2766 


} 



(i) 

Table 1: Contributions to the pressure at 9 = for Bohm backgrounds. P x = kiH c 



5 Summary, discussion and open questions 

Large generalised black holes with a negative cosmological constant tend to have dynamical 
instabilities if the horizon has a curvature scale which in some region is large compared to 
the constant Ricci scalar of the horizon, V\ 3> R. This predicts a critical temperature in 
the dual field theory: Tq ~ OK. 

We specialised to the case of a six dimensional field theory, as this corresponds to the 
lowest dimension in which unstable black holes are known and further may be considered 
within the AdSj x S 4 version of the AdS / CFT duality. 

We have argued that field theories on curved backgrounds with a curvature scale OK S> R 
can generically have a thermodynamic instability at Tq ~ 9^ due to the strong curvature 
region having negative compressibility. The instability is caused by contributions to the 
pressure from mass dimension six curvature tensors. These come from mass dimension 
curvature six scalars that are renormalisable vacuum terms in the field theory action. We 
calculated the explicit form of the contribution to the energy momentum tensor, and eval- 
uated these tensors for the case of Bohm backgrounds. 

There are various unsatisfactory aspects of our treatment. We do not have a precise 
matching. Whilst a precise criterion for black hole instability is known [22], this is not 
the case in field theory. This is because the coefficients of the curvature terms in the 
energy momentum tensor are not predicted by the field theory, but are rather empirically 
determined, just like coupling constants. Yet, the instability itself depends on the total 
contribution to the pressure having a specific sign. 
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This situation potentially raises a conceptual problem for the AdS/CFT correspondence. 
The field theory has undetermined coefficients that correspond to finite terms in the bound- 
ary action on the gravitational side of the duality. However, these terms do not influence 
the dynamics of the bulk, and the classical stability of the bulk in particular. A possible 
resolution of the problem is that the bulk spacetime determines a specific renormalisation 
scheme for the field theory terms in question. In order for the matching of instabilities to 
be made precise, this renormalisation scheme would have to be known explicitly. 

Nonetheless, the qualitative agreement is suggestive. It fits well into the recent discovery 
of agreement between thermodynamic and classical instabilities of horizons, and extends 
the connection to cases without translational invariance. An exciting possibility is that 
the ideas presented here could be adapted for use in other contexts where the stability of 
inhomogeneous horizons is of interest. In more general cases the Casimir energy would 
presumably not be that of a dual field theory, but rather due to gravitational degrees of 
freedom. A possible application of current interest is the stability of nonuniform black 
branes [8, 9, 10, 11]. Can the stability be studied from a thermodynamic perspective? Can 
one associate a Casimir energy and pressure with general inhomogeneous horizons? 
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A Divergences of mass dimension six two-tensors 

This appendix contains formulae that are necessary for finding the independent linear com- 
binations of the tensors in (23) that have vanishing divergence. This is done by expressing 
the divergences of each rank two tensor in terms of independent rank one tensors. The 
independent rank one tensors will be 

AP = V a R pqrs V u V a R™ rs , 

a(2) _ ft n q ^janpers 
— ±L pqrs ±L e av v 11 > 

4(3) _ f> B 1 r X7 E epas 
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4(4) — U D q yanpers 



(26) 



With a little work, one may then show that 





► +0( J R9i 5/2 ), 



(27) 




1 .(2) _ 1.(3) 

-A^ + 2A^ 



where the 0{R^ > / 2 ) indicates that throughout the calculation, all Ricci tensors have been 



dropped. The results are exact if R^ v = 0. Recall that we do this because we are interested 
in the dominant contributions in the region where !H> iJ. 

There is a fifth tensor that appears in these calculations which one might think is linearly 
independent 



the backgrounds we are considering, R pqrs R pqtu R rs tu is not 0(9i 3 ). We are only keeping 
track of leading order terms, so this term is not linearly independent. 

Given these expressions for the divergences in terms of linearly independent one-tensors 
(27), it is simple linear algebra to find the two divergence free tensors of equation (24). 
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